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ABSTRACT: In this paper we study qualitative features of glueballs in M = 1 SYM for 
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spectrum without any cut-off. We propose a regularization procedure needed to make these 
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1 Introduction and General Ideas 



During the sixties the study of hadronic particles was a mainstream area of theoretical 
Physics. The Regge trajectories were proposed and by the systematic study of scattering of 
Hadrons, the dualities between the s and t channels amplitudes were investigated and nicely 
realized by the Veneziano amplitude, 

The dual models, precursors of the modern String Theory, were one of the best candidates 
to explain the relevant Physics and provided the tools to explore the issues mentioned above. 
Despite some success, the experiments realized during the late sixties involving scatterings 
of large s and large t Mandelstam variables, but keeping fixed | (fixed angle processes), gave 
results that showed that the amplitude was falling as a power law instead of the exponential 
law predicted by the dual models. This situation, together with the advent of QCD (that 
correctly predicted the scaling for fixed angle scattering and many other things) , lead to the 
demise of the dual models for the study of hadronic Physics. 

Even when the original motivation to study dual models momentarily dissapeared, their 
rich structure kept many physicists interested and, after some technical subtleties were un- 
derstood, this was the beginning of String Theory. 

As is known, after thirty years, string theorists have come back to the study of problems 
related to the hadronic world. Indeed, guided by the Maldacena Conjecture pQ and their 
refinements j2], there have been very many interesting achievements in the area. They 
correspond to field theories with different amount of SUSY (including no SUSY) but "very 
similar" to QCD. The important point is that many characteristic features of QCD, like 
confinement, chiral symmetry breaking, etc; have been understood based on dual String 
theory backgrounds. 

In this paper, we study glueballs in some of the models mentioned above as being "very 
similar" to QCD (even when the models we will deal with here and those available in the 
literature, perhaps are not in the same universality class of QCD). Let us motivate a little 
bit the study of these glue-composed excitations. 

We know that the main distinction between a field theory in a confining phase and the 
same field theory in the Higgs phase is the presence of Regge trajectories, that do not occur 
in theories with Coulomb of Yukawa interactions. These Regge trajectories appear when 
plotting the spin J and the squarred mass m 2 of the excitation, thus giving relation of 
the form J = a'm 2 + a Q , with a' ~ (lGeV)~ 2 ; this relation does not have in principle, 
an upper bound in J, m. It is due to these infinite number of Regge resonances, being 
interchanged in the s, t channels of any hadron scattering that the beautiful structure of 
duality appeared in the models above mentioned. The glueballs should be some of these 
Regge excitations (making up a full trajectory if mixing with quarks is neglected) and this 
is a possible motivation to study them. 

From a modern QCD perspective, it is known that the cloud of gluons is what logically 
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connects between a current quark (with mass of a few MeV) to a constituent quark, with 
mass of around 300 MeV. Since glue is part of the hadronic matter, we can consider color 
singlet composites of the form qqg, gg, ggg (apart from the mesons, baryons and exotics). 
The glueballs are composites made out of constituent glue, with no quark content. Of course, 
since we live in a world with quarks, one might think that the proposal of pure glue objects 
is impossible to study, because the quarks should run in loops when doing corrections to the 
operators, that leads to glueballs mixing with mesons, rendering the object not-pure glue. 
But lattice theorist (working in the quenched approximation) are not stopped by this. Indeed 
they took advantage of the limitation and have taught us many things about glueballs. 

Among the things that Lattice showed about QCD glueballs, we can mention the facts 
that: 

• there is a bound state spectrum 

• the lightest glueball is a scalar 

• the next is a tensor, 1.6 times heavier 

• the mass of the lightest glueball should be around 1630 MeV; 

see for example j3] for a nice and clear review of these results. 

How are these lattice predictions experimentally checked? Experimentalist look for pro- 
cesses rich in glue production, like the J / if) decay, where the c, c quarks annihilate into 
gluons. Other process might be the pp annihilation, in this case the idea is that the quarks 
and anti-quarks in the initial hadrons annihilate completely, producing glue that later decays 
into hadrons. There are many glueballs candidates. One of them seems to be well established 
and is called / (1500) with a width of 112 MeV g]. 

From the string theory view point, using the Maldacena duality for the case of confining 
backgrounds, the study of glueballs (in the field theory dual to the background) proceeds 
by finding bound states for the fluctuations of the supergravity fields. Basically, the idea 
is to fluctuate all the fields in a given IIA, or IIB solution dual to a confining field theory, 
and linearizing in the fluctuated fields, study their eqs of motion (that are the Einstein, 
Maxwell and Bianchi eqs). The system is reduced to a Schroedinger problem. When solved, 
has eigenf unctions that we identify with the glueballs and eigenvalues that are identified 
with their masses. The fluctuations of the fields are dual to different operators in the gauge 
theory and what we are actually computing in the field theory side is the two-point correlation 
function of two glueball operators that should behave in a Wilson expansion as 

< 0{x)0{y) - '' " 

j 

where Mj are the glueball masses. The quantum numbers of the glueballs J PC are determined 
on the basis of the spin (J) of the supergravity field and the R-symmetry quantum numbers 
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(in a KK-harmonics decomposition) as studied, for example, in j^j. We should point that 
this procedure is not totally clear in many of the available confining-models and it should 
be important to understand it better. 

This machinery has been applied to some confining models. Let us add that, since many 
of the existing Supergravity models are duals to confining field theories with only adjoint 
matter content, the objects under study are only glueballs (no hybrids) and since we work 
in the large N c regime, the glueballs are stable. 

Let us briefly review what was done in this subject, the original idea, described above, has 
been proposed by Witten in [Hj. Many papers followed, exploring this nice idea in different 
contexts. For example, confining models using black hole geometries were developed for 
QCD 3 and QCD4 in [Zj, [Hj. Also, models based on rotating branes were introduced [9 j and 
other models based on AdS§ with no SUSY ^0]- All these models have an spectrum that is 
numerically very close to the one obtained by Lattice methods. 

We should stress, that even when the comparisons between the lattice and "AdS" based 
results seem so accurate and promising, these calculations are done in opposite regimes. 
Indeed, the gravity-dual computation is in strong 't Hooft coupling and this limitation is 
imposed by the Supergravity approximation. On the other hand, the Lattice computations 
are done at weak coupling, this seems to be a necessity of having a continuum limit because 
the lattice spacing a, has a relation clKqcd ~ e~ x l9 lN with the QCD coupling and scale 
(other regularizations give similar results). One might think about doing strong coupling 
lattice computations, but they do not seem to be smoothly related to the continuum theory. 
It is possible that the numerical coincidences aluded above, are based on some dynamical 
principle to be understood. 

There exist a set of Supergravity models that preserve N = 1 SUSY that have been object 
of lots of study and amusing advances. One of the models was put forward by Klebanov 
and Strassler jTTj and the glueballs in this model were carefully studied in the set of papers 
|12j.|13|. The results indicate that, for the Klebanov- Strassler model, the masses of the 
++ , 1 , 2 ++ in the strong 't Hooft coupling limit, fall in a linear trajectory. Also, the 
finding of a massless excitation showed that this cascading field theory is not in the same 
universality class of QCD, because of the reasons we explained above. 1 

1.1 Motivations and organization of this paper 

We mentioned above a set of Supergravity duals to confining models and up to this point 
we just commented on the one proposed by Klebanov and Strassler [TJ]. There exist some 
other models that are based on wrapped D-branes. The main idea here is to consider the low 
energy field theory in (k + 1) dimensions, obtained by wrapping a Dp brane on a {p — k) cycle. 
Some subtleties of this type of models will be explained in section 2 of this paper. Here we 

■""We thank Oliver Jahn for extensive discussions on many of the points touched in this introduction 



3 



just want to emphasize that the glueballs spectrum in this case is poorly understood. 

Indeed, there is a paper [H] , where a study was initiated. We believe that this study is not 
completely correct from a technical viewpoint (we believe that incorrect eqs were used) and 
the conclusions expressed there are, even when intuitively understandable, also not totally 
correct. The main point of that paper is that in one of these models it is necessary to 
introduce a hard cut-off in order to have a discrete spectrum. In this paper we re-analize 
this statement and propose a different result, basically that all these models do have discrete 
spectrum of glueballs and give a way of computing it. The spectrum even though discrete is 
not normalizable. In order to get normalizable states we need to introduce a regularization 
procedure. The regularization we propose is not the introduction of a hard cut-off but is 
more in the spirit of the Wilson loop calculations where a non-physical part is 

subtracted. In this paper we will not make much emphasis on the numerical aspects of the 
problem. Indeed, even when discrete states are numerically obtained, we will not worry 
here about comparisons with lattice results, that as explained above are perhaps not very 
significative. The main objective of this work is to study qualitative features of the spectrum, 
point out differences with previously studied cases, and propose a procedure of computing 
and regularizing in these wrapped branes set-ups. 

This paper is divided in two parts, one dealing with a particular type IIB model and the 
other with a type IIA model. Both parts have been written and can be read in parallel and 
almost independently. 

In section 2, we describe in detail the two models we will be using, one based in type 
IIB, with D5 branes wrapping a two-cycle inside the resolved conifold. The other in type 
IIA, based on D6 branes wrapping a three cycle in the deformed conifold. Section 3 deals 
with the glueballs in the type IIB set-up, while section 4 sketches the results corresponding 
to the type IIA model. We did not carry the problem to an end because of the need of 
a more precise numerical analysis, since the solution is only numerically known. Section 5 
presents conclusions and possible future work proposed to the interested reader. There are 
very detailed Appendixes that carefully explain all the computations in sections 3 and 4. 

2 J\f = 1 SYM models from wrapped branes 

In this section, we write an account of duals to N=l SYM from wrapped branes. The two 
main models on which we will concentrate are the ones based on D5 branes wrapping a two 
cycle, that we will consider to be a two-sphere inside a CY3 fold and D6 branes on a three 
cycle (a squashed three sphere), also inside a CY3 fold. We will present the solutions in 
detail, and explain the main characteristics of the dual gauge theory. We will emphasize the 
existence of 'extra' modes called KK modes with mass of the same order of the confinement 
scale. Since our interest in this paper is on glueballs, we will discuss the influence of this 
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'extra' modes in the computation of glueballs for N = 1 SYM. 



2.1 D5 branes wrapping S 2 

We will work with the model presented in ^H] (the solution was first found in a 4d context 
in [TH]) and described and studied in more detail in the paper [20] • Let us briefly describe 
the main points of this supergravity dual to AT = 1 SYM and its UV completion. 

Suppose that we start with N D5 branes, the field theory living on them is (5 + 1)SYM 
with 16 supercharges. Then, suppose that we wrap two directions of the D5 branes on a 
curved two manifold that can be choosen to be a sphere. In order to preserve SUSY a 
twisting procedure has to be implemented. The one we will be interested in this section, 
deals with a twisting that preserves four supercharges. In this case the two-cycle mentioned 
above lives inside a CY3 fold. Notice that this supergravity solution will be dual to a four 
dimensional field theory, only for low energies (small values of the radial coordinate). Indeed, 
at high energies, the modes of the gauge theory start to explore the two cycle and the theory 
becomes first N=l SYM in six dimensions and then, the blowing-up of the dilaton forces us 
to S-dualize and a little string theory completes the model in the UV. In this sense, to study 
only the 4d-SYM part of the background, a procedure that "substracts" the unwanted UV 
completion, should be useful. We will elaborate on this in Section 3.1. 

The supergravity solution corresponding to the case of interest in this section, the one 
preserving four supercharges, has the topology of R 13 x R x S 2 x S 3 and there is a fibration 
between the two spheres that allows the SUSY preservation. The topology of the metric, 
near r = is R ' 6 x S 3 . The full solution and Killing spinors are written in detail in [20J. 
Let us revise it here for reference. The metric in Einstein frame reads, 

ds 2 = a'g s Nei \ —^—dx\ 3 + e 2h ( d9 2 + sin 2 Od<p 2 ) + dr 2 + \ U - A') 2 1 , (2.1) 
l a'g s N ' A 

where <fi is the dilaton. The angles 9 G [0, tt] and ip e [0, 2ir) parametrize a two-sphere. 
This sphere is fibered in the ten dimensional metric by the one-forms A 1 {% — 1, 2, 3). Their 
expression can be written in terms of a function a(r) and the angles (9, ip) as follows: 

A 1 = -a(r)d6, A 2 = a(r)sm9d(p, A 3 = - cos 9d<p. (2.2) 

The w % 's appearing in eq. (|2.1|) are the su(2) left-invariant one-forms, satisfying 

w 1 = cos ipd9 + sin ip sin 9 dip , 

w 2 = — sin ipd9 + cos ip sin 9d(p , 

w 3 = dtp + cos9d(p . (2.3) 
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The three angles (p, 9 and ip take values in the rank < ip < 2tt, < 9 < it and < ip < An. 
For a metric ansatz such as the one written in (|2.1j) one obtains a supersymmetric solution 
when the functions a(r), h(r) and the dilaton </> are: 



air) = — , 

smh 2r 

e 2h = rcoth2r - 



-24> _ o -20 o 



sinh 2 2r 4 ' 

2e h 



sinh2r ' ^ 

where 0o is the value of the dilaton at r = 0. Near the origin r = the function e 2h behaves 
as e 2h ~ r 2 and the metric is non-singular. The solution of the type IIB supergravity includes 
a Ramond-Ramond three-form given by 

^F (3) = - 1 -(w 1 -A')A(w 2 -A 2 )A(w :i -A 3 ) + ~ ^ F a A (w a — A a ) , (2.5) 

where F a is the field strength of the su(2) gauge field A a , defined as: 

F a = dA a + ^e abc A b AA c . (2.6) 
The different components of F a are: 



F 1 



-a' dr A d9 



a' sin 9dr A dip 



1- a 2 ) sm9d9 A dip , (2.7) 



where the prime denotes derivative with respect to r. Since dF^ = 0, one can represent F( 3 ) 
in terms of a two- form potential C(2) as Frg\ = dC^)- Actually, it is not difficult to verify 
that C(2) can be taken as: 



(2) 



N 



1 r 



ip ( sin 9d9 A dip — sin 9d9 A dip ) — cos 9 cos 9dip A dip — 



-a(d9 Aw 1 - sin 9dip A w 2 ) 



(2.1 



Moreover, the equation of motion of F^ in the Einstein frame is d(^e^*F^ ) = 0, where * 
denotes Hodge duality. Let us stress here that the previous configuration is non-singular. 

Finally, let us comment on the fact that the BPS equations also admit a solution in 
which the function a(r) vanishes, i.e. in which the one-form A 1 has only one non- vanishing 
component, namely A 3 . We will refer to this solution as the "abelian" (or "singular") M = 1 
background. Its explicit form can be easily obtained by taking the r — > oo limit of the 
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functions given in eq. (|2.4|) . Notice that, indeed a(r) — > as r — > oo in eq. (|2.4p . Neglecting 
exponentially suppressed terms, one gets: 




while can be obtained from the last equation in (|2.4jl . The metric of the abelian background 
is singular at r = 1/4 (the position of the singularity can be moved to r = by a redefinition 
of the radial coordinate). This IR singularity of the abelian background is removed in the 
non-abelian metric by switching on the A 1 , A 2 components of the one-form (|2.2jl . 



2.2 Some analysis of this model 

Let us first summarize the field theory aspects of the dual to the gravity solution we will 
be mainly concerned with. The main characteristic is that it contains a four dimensional 
Minkowski space, a radial direction and a two sphere fibered over a three sphere. In j^l this 
solution was argued to be dual to N — 1 SYM. Let us analize the claim a little more, the 
field theory at low energies (low compared to the inverse size of the two-sphere) has degrees 
of freedom given by a vector field and a Majorana spinor (in 4d). When increasing in energy, 
other modes with mass of the order of the inverse size of the 5 2 appear in the spectrum. 
These are called KK modes and can be seen as coming from the reduction of the maximally 
D5 branes SUSY field theory on a two dimensional sphere and a twisting (explained below) 
are performed. When the energy is high enough the excitations of the theory propagate 
in 5 + 1 dimensions and the UV-completion of our minimally SUSY four dimensional field 
theory is the six dimensional little string theory living on N NS5 branes. 

Let us recall briefly the twisting procedure. One has a 5 + 1 field theory (that lives on 
N D5 branes) that has gauge fields, fermions and four scalars, all in the adjoint of the 
SU(N) gauge group. We rewrite the 50(1,5) x SU(2)l x SU(2)r group quantum numbers 
of the fields above in terms of 50(1,3) x 50(2) x SU{2) L x SU{2) R and then we mix the 
quantum numbers respect to 50(2) with those of another 50(2) that lives inside one of 
the SU(2)'s. After this twisting procedure is performed, we are left with fields that under 
50(1,3) x U(l) x SU(2) transform as 

^ = (4,0,1), $ a = (l,±,l), a = 2(l,±,2). (2.10) 

for the bosons that can be seen to be a massless gauge field, a massive scalar (coming from 
the gauge field) and other massive scalars (that originally represented the positions of the 
D5 branes in R 4 ). As a general rule, all the fields that do transform under the twisted U(l), 
the second entry in the charges above, will be massive. For the fermions we will have 

r = (2, 0, 1), (2, 0, 1), (2, ++, 1), (2, — , 1), (2, 0, 2), (2, 0, 2) (2.11) 
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that is a Majorana spinor in four dimensions that is massless and then we have massive ones 
(those whose quantum number under the twisted U(l) is not zero). The KK modes are the 
massive modes mentioned above. Their mass is of the order M\ K = (Rs2)~ 2 = g ^, N - 

The dynamics of these KK modes, mixes with the dynamics of confinement in this model, 
because the scale of strong coupling of the theory is of the order of the KK mass. If we could 
work with a sigma model for the string in this background (or in the S-dual NS5 background) 
to all orders in a' we could decouple both behaviours. Meanwhile, the dynamics of these KK 
modes has not been studied in great detail, but some progress have been made, for example 
in the papers [21]. Finally, we would like to mention a paper where a very careful study of 
the KK modes spectrum have been done, also pointing a coincidence with N = 1* theory in 
a given Higgs vacuum [22J. 

Let us briefly comment on the influence of these KK modes on the glueballs spectrum. 
Indeed, once the strong coupling regime of the field theory is attained, one possible way to 
compute is using these supergravity backgrounds. Given that we are in the supergravity 
approximation, the spectrum of our model includes these KK 'contaminations' (this feature 
repeats in all the dual to non-conformal field theories). Obviously, our glueballs will be of 
two types, those coming from condensates of the gluon and gluino, those 'composed' out of 
KK modes and finally, hybrids, composed out of SYM fields and KK modes. We would like 
to discard those with some KK constituent. We will comment in the conclusion section on 
a possibility to do this. 

Finally, let us mention that there are many succesful checks showing that the supergravity 
background presented above captures different non-perturbative aspects of N — 1 SYM. We 
will not discuss these many checks here, instead, we refer the interested reader to the very 
careful reviews [23J. 

2.3 D6 branes wrapping 

Now, let us comment on the models based on D6 branes wrapping a calibrated three-cycle 
inside a CY3 fold. The progress in this direction originated from the duality between Chern- 
Simons gauge theory on S 3 at large N and topological string theory on a blown up Calabi-Yau 
conifold 25]. This duality was embedded in string theory as a duality between the IIA string 
theory of N D6-branes wrapping the blown up S 3 of the deformed conifold and IIA string 
theory on the small resolution of the conifold with iV units of two form Ramond-Ramond 
flux through the blown up S 2 and no branes [2Ej. The D6-brane side of the duality involves 
an M = 1 gauge theory in four dimensions that is living on the non-compact directions of 
the branes, at energies that do not probe the wrapped S 3 . 

Just like before, in order for the wrapped branes to preserve some supersymmetry, one 
has to embedd the spin connection of the wrapped cycle into the gauge connection, which is 
known as twisting the theory. 
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When we have flat D6 branes, the symmetry group of the configuration is 50(1,6) x 
SO(3)r. The spinors transform in the (8,2) of the isometry group and the scalars in the 
(1,3), whilst the gauge particles are in the (7,1) j2Zj- Wrapping the D6 brane on the three- 
sphere breaks the group to £0(1,3) x £0(3) x £0(3)^. The technical meaning of twisting 
is that the two £0(3)s get mixed to allow the existence of four dimensional spinors that 
transform as scalars under the new twisted £0(3) [2H]- One can then see that the remaining 
particles in the spectrum that transform as scalars under the twisted £0(3) are the gauge 
field and four of the initial sixteen spinors. Thus the massless field content is that of M = 1 
SYM. Like in the model analyzed in the previous section, apart from these fields, there 
will be massive modes, whose mass scale is set by the size of the curved cycle. When we 
probe the system with very low energies, we find only the spectrum of M = 1 SYM. For D6 
branes in flat space, the 'decoupling' limit does not completely decouple the gauge theory 
modes from bulk modes In our case, we expect a good gauge theory description only 
when the size of the wrapped three-cycle is large, which implies that we have to probe the 
system with very low energies to get 3+1 dimensional SYM (3U]- In this case, the size of the 
two cycle in the flopped geometry is very near to zero, so a good gravity description is not 
expected. In short, we must keep in mind that the field theory we will be dealing with has 
more degrees of freedom than pure M = 1 SYM, thus the glueballs masses that one might 
obtain following the procedure explained in the following sections might be 'contamined' by 
glueballs composed out of KK modes or hybrids composed out of KK modes and gluons 
or gauginos. Again, how to decouple the ones we are interested into from those glueballs 
'composed' of KK modes is going to be discussed in the conclusion section. 

Finally, let us add that the duality described above is naturally understood by considering 
M-theory on a G2 holonomy metric [SHj. In eleven dimensions, G2 holonomy implements 
M = 1 as pure gravity. One starts with a singular G 2 manifold that on dimensional reduction 
to IIA string theory corresponds to iV D6 branes wrapping the £ 3 of the deformed conifold. 
There is an SU (N) gauge theory at the singular locus/D6 brane. This configuration describes 
the UV of the gauge theory. As the coupling runs to the IR, a blown up £ 3 in the G 2 manifold 
shrinks and another has fixed size. This flop is smooth in M-theory physics. The metrics 
will be discussed in more detail in the following sections. In the IR regime, the G2 manifold 
is non-singular and dimensional reduction to IIA gives precisely the aforementioned small 
resolution of the conifold with no branes and RR flux. 

Let us now, write explicitly the background on which we will be interested. It is conveninet 
to start with the eleven dimensional M-theory background, that reads 

ds 2 n = dxl 3 + ds 2 7 (2.12) 

with 

ds 2 = dr 2 + a(r) 2 [ (Ei + g{r)ai) 2 + (Ea + g (r )(T2 y-\ + c(r) 2 (S3 + g 3 ( r ) a3 f 

+&(r) 2 [a 2 + a 2 2 l+/(r)V 2 , (2.13) 
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where Ej,(7j are left-invariant one-forms on the SU(2)s f)2.3j) . The six functions are not all 
independent 

. , —a(r)f(r) , . , , 9 , 

g(0= 26(r)c(r) » ^(r) = -l + 2^(r) 2 . (2.14) 

None of the radial functions are known explicitly, although the asymptotics at the origin and 
at infinity are known. The asymptotics are found by finding Taylor series solutions to the 
first order equations for the radial functions. The equations are PS].|3§] 

a ~ 2a + 86 4 c 3 ' ~ 2b 8b 3 c 3 

c 2 c 2 _ 3a 2 f 2 ■ _ _ a 4 / 3 

C "~ 1 + 2^ + 2p"^6 4_ ' /__ 4W (2 - 15) 
As r — > one has 

r _ (g 2 + 2)r 3 _ (-74 - 29g 2 + 31g 4 )r 5 
^ ~ 2 288i?g 69120.R 4 ' 

, o (9o 2 ~ 2)r 2 (13 - 21q 2 + llg 4 )r 4 

b{r) ~ R °-^QR TT52RI + ■■■' 

r _ (5g 2 - 8)r 3 _ (232 - 353g 2 + 157g 4 )r 5 
~ ~2 288RI 3456UR^ + '"' 

where go an d -Ro are constants. Note that a(r) and c(r) collapse and the other two functions 
do not. As r — > oo we have 

r v^gi^i (27^6 - 96^i)i?? 

a(r) = —pz = 1 1 , 

1 ; v/6 V2 96r 

Kr) _ + ¥ + _ 

-r „ 9i? 2 

With constants Ri,qi,hi. Note that /(r) stabilises. Three constants appear to this order, 
whilst there were only two constants in the expansion around the origin. This just means 
that for some values of these constants, the corresponding solution will diverge before it 
reaches zero. In any case, we find no hi dependence in the results below. 

We can reduce this to Type Ha and we will find a non-singular background with dilaton, 
metric and RR one form excited, that reads, 

ds 2 IIA>stHng = 2e 2 ^(dxl 3 + dr 2 + b{rf{a{ + a 2 2 ) + a(r) 2 ((£ x + g^f + (E 2 + g(r)a 2 ) 2 ) + 
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f 2 c 2 



4e 4/3^ = /(r) 2 + c(r)2(l + (r))2 ^ = cog ^ + cog ^~ + / I ^ 3 - E 3 )(2.18) 

f 2 + c 2 {l + g 3 ) 2 

Where we have defined dxn = dip' + dip and dip = dip — dip'. So, to summarize the things 
clearly, let us write the metric of our IIA solution in Einstein frame (as will be used below), 

ds% = 2e* /6 (dxl 3 + dr 2 + (b(r) 2 + g(r) 2 )(d9 2 + sin 2 9d<p 2 ) + a(r) 2 (d9 2 + sin 2 9d(p 2 ) + 
+2 g(r)a(r) 2 [cos ip(d9d9 + sin 9 sin 9d<pd<p) + simp (sin 9d9dip — sin9d9d0)] + 

x 2 2 

+ « {?7 ^(# + cos^-cos^) 2 ) (2.19) 

/ 2 + c 2 (l+^ 3 ) 2 ' 

with the same dilaton as in (|2.18j) . besides, the field strength F 2 reads, 

F 2 = k'(r)drA(dip+cos9dif-cos9d0)-(k(r)+l) sm9d9Adp+(k(r)-l) sm9d9Ad<p, (2.20) 
and 

/^-^(i- g j) 

U / 2 + c2(l+ 93 P' 

To end this section, let us briefly revise what checks exist of the duality between the 
backgrounds presented here and N — 1 SYM. 

Of course, the number of supercharges match, there is a nice picture of confinement in 
terms of a Wilson loop computation in IIA. But most of the presently known matchings 
of N = 1 SYM with G 2 holonomy M-theory come from considering membrane instantons 
as gauge theory instantons that generate the superpotential (HI], membranes wrapped on 
one-cycles in the IR geometry that are super QCD strings in the gauge theory 155] . 
and fivebranes wrapped on three-cycles that give domain walls in the gauge theory |32| I34j . 
These matchings above, are essentially topological and do not use the explicit form of the 
G 2 metrics. In the category of test/checks that use the form of the metric, we can mention 
|35j . where rotating membranes in these G 2 backgrounds have been studied and relations 
for large operators in SYM have been reproduced. We should also mention jHE], where a 
very nice picture of the chiral anomaly of SYM have been developed. Perhaps less promising 
is the fact that confining string tensions do not arise as cleanly in the IIA backgrounds as 
in the type IIB case [HI!- There are many aspects of this duality that are not on a very 
firm basis and we think that through study, these unclear points might become clear. The 
results that we will present in Section 4, use the explicit form of the metric and should be 
considered to belong to this second category of tests. 



11 



3 Glueballs from type IIB solution 



As we explained above, in order to study glueballs, we need the variation of the eqs of motion. 
In the type IIB case, for the solution of D5 branes wrapping S 2 inside a Cl^-fold discussed 
in the previous section, the Einstein eqs for the metric, dilaton and three form read, 2 

irW = \d^ v <j> + 9 -^-(F^F* v * - ^g^Ff) (3.21) 
by contracting we get the Ricci scalar eq, 

R = \{d»<P?+ 9 -^-Fl (3.22) 
The dilaton, Maxwell and Bianchi eqs. are, 

V 2 = 9 -^-Fl d^y/getF*"") = 0, d [lx F Rvp] = (3.23) 

Now, let us study the fluctuations of these eqs. Let us assume that the background fields 
vary according to 

9^ -»• g^ + eh,,,,, </> -> + e5(p, F^p F lu/p + eSF^p (3.24) 

Keeping only linear order in the parameter e, we get eqs for the fluctuated fields (that can 
be found written in detail in Appendix B). The metric fluctuation, can be splitted in its 
worldvolume, internal and mixed parts: h^, h a $ and hi a respectively. We assume that there 
is no fluctuation in the mixed part, h ia = and perform a Weyl shift [S] in the worldvolume 
fluctuation, 3 

hi = Kj + -^gijKi » Kp = h (af3) + \g a pK (3.25) 

We denote with latin indices the worldvolume and transverse coordinates i,j = x,t,r and 
with with greek indices the internal coordinates, a, /3 = 9,9,(p,{p,ip. Also, A is a constant 
and h( a p) is the symmetric traceless part of the fluctuation in the internal directions. Notice 
that the trace of the metric fluctuation over all the space is h)t = (A + From here 

on we will denote with x the radial and 'gauge theory' coordinates and with y the internal 
coordinates (the angles on the spheres). Imposing a de Donder and Lorentz type condition, 
V a /i( a/ 3) = V a h ai = 0, the decomposition in harmonics for the fluctuated fields is, 

K{x,y) = J: I U I (x)Y I (y), h (aP) = E I b I (x)Y^ ) (y), 

hij = EjH l3 {x)Y\y), 5<j>{x,y) = i: I f I (x)Y I (y) 

Cjj = S ja f 3 -(x)y f (y), Cjg = S f qf(x)y Q J (y), C a/3 = E I a I (x)Y^(x). (3.26) 

2 We will denote the contraction of indexes with * symbols, so, for example g^ v A^iB^jp — ApkiBjp = 
A*kiB* p . 

3 The standard Weyl shift in a D dimensional spacetime is A/5 = —1/(1? — 2); this value is required to 
simplify the variation of V 2 . Here we choose to leave A as a constant to be determined later. 
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Where C2 is the two form potential. We set h'- = 0. Given this ansatz the equations 
of motion can be consistently solved for the other fluctuations. Using the decomposition 
in harmonics, keeping only the s-wave and choosing a particular value for A = — 1| the two 
equations for the fluctuations nicely combine into just one equation (again, see the Appendix 
B for details). This final eq. reads, 



9» 
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29-3, ab „ 35 



V 2 f(x) + 2g rr d r f(x)d r cj ) + ^e* -Fi + ——(g" b F a ^F b ** - -g rr F m F r **) f(x) = (3.27) 



Expanding (|3.27|l in plane waves, f(x) = F(r)e^ K ' x \ we have, 



dr 2 \ yfg J dr 

+ + ^(tfFv-FT - %9 rr F r **F?)\ - KY X ) F(r) 



(3.28) 



This equation will be solved numerically; the glueball masses are given by the eigenvalues 
K 2 for which there is a solution with appropriate boundary conditions. Before studying the 
boundary conditions of this problem let us cast eq. (J3.28|) in a more familiar way. To save 
us some writing, denote the coefficient of the first derivative term (^) in ()3.28|) as fj,(r) and 
the coefficient of F(r) as a(r), that is, 

Ai(r) = 2d^ + ^Ld r {^) 

a{r) = 9 ^e*fy<$ + 2 -^(9 ab F^Fr - (3.29) 

Making a change of variables, J~(r) = e~zf ^ dz F(r), equation ()3.28j) can be written in a 
Schrodinger form 

d?F(r) 
dr 2 

where 



VS{r)F = 0, (3.30) 



VS(r) = -a{r) + i(/i 2 (r) + 2^^) + K 2 (3.31) 

A graph of the Schrodinger potential VS(r) is given in Figure [0 At this point it is convenient 
to recall that, as explained in section two, the full model, considering its UV completion is 
not dual to a non-abelian gauge theory, but to a little string theory. Indeed, in the UV, 
due to the divergent dilaton, the solution has to be S-dualized yielding a IIB solution with 
NS five branes. In the decoupling limit this background is not dual to M = 1 SYM in 
four dimensions but to a higher dimensional, 5+1, little string theory. Naturally, trying to 
calculate observables by simply using the solution up to infinity might not yield sensible 
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Figure 1: VS(r) for different values of K 2 . 



answers. In what follows we want to study the glueball spectrum and, if necessary, propose 
a regularization procedure. 

To calculate the ++ mass we have to numerically find eigenvalues satisfying equation and 
appropriate boundary conditions. The asymptotic behavior of the potential is, 

Q2 

VS(r -> oo) = -— + K 2 (3.32) 
5 

VS (r^0) = -^ + K 2 (3.33) 
Choosing the exponentially decreasing solution at infinity we get, 



T{r -> oo) ~ e -V-f+ K2r . (3.34) 

At the origin we demand a smooth solution, 

dF(r) 



dr 



0. (3.35) 



Similarly to Klebanov-Strassler, satisfying the boundary conditions implies that the eigenval- 
ues are bounded from below, K 2 > 92/5 and thus, there is a mass gap. But here, in addition 
of being bounded from below, the eigenvalues are also bounded from above, K 2 < 1036/45. A 
similar phenomenon was observed in |30J in the context of non-commutative gauge theories. 

Another difference with the Klebanov-Strassler model is that here the boundary condition 
at infinity depends on the eigenvalue. In a technical sense, each eigenvalue defines a different 
problem -different boundary condition- and the spectrum is then given by the eigenvalues 
of this collection of problems; It is a more general situation than the standard eigenvalue 
problem. 
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Using the WKB method ^2] we can estimate the eigenvalues K 2 for which there exists 
a solution of ()3.30|) satisfying ()3.34j) and ()3.35j) . Also, it can be shown numerically that 
the WKB integral f^ K ■* \J—VS(r)dr is a monotonically decreasing function of K 2 . And 
this fact can be used to prove that there is only one eigenvalue in the spectrum. We find, 
K = 4.33 (1/ y/g s a'N). However, it is easy to show that this eigenvalue does not correspond 
to a normalizable state. For large r, ^det g ~ e^ r and thus, 

y fi^jjfi „ c (|-V(-¥+^>) (3.36) 

so the integral / °° ^detgpQ does not converge. The issue we are confronted with now is to 
find a good regularization for this model. 

Let us note an important point. Sean Hartnoll pointed that one might think of taking 
the norm in flat space / c?r|jF| 2 . Indeed, the equation ()3.30|) seems to indicate that, but we 
used a norm obtained from on the ten dimensional curved background / <i 10 x|jF| 2 , and is 
this norm the one that is forcing us to some regularization. His comment is based on the 
fact that one should only worry about our fluctuations to have finite energy and according 
to the paper |Hj, this condition implies the finitness of the norm in flat space. If we use 
his proposal, our wave functions T{r) are normalizable. In this paper we choose to attach 
to the more conventional norm defined in the curved space. The differences and physical 
implications of each choice of norm will be investigated elsewhere. 4 

Coming back to the regularization we need to introduce, let us recall the reader that, as 
explained in section two, in addition to the solution given by (|2.1|) - (|2.8j) there is another 
solution to the BPS equations, that we aluded to as the "abelian" solution. The abelian 
(singular) and non-abelian (non-singular) solutions are the same at infinity but the abelian 
solution is singular in the IR and thus is not a good dual to SYM, that has nothing singular. 
We propose that a good regularization for this model is to use the background presented in 
(I2.1j) - (j2.8j) up to the point where it becomes indistinguishable from the abelian one. After 
this point the two solutions are the same and neither of them is dual to SYM, but to a higher 
dimensional theory on which we are not interested here. 

Figure (j2J shows a plot of the Schrodinger potential for the abelian and non abelian 
solutions. Our proposal is that the IR solution that captures the physics of N=l SYM is 
valid only up to the region A^b where the potential becomes the same as the one of the abelian 
solution. The scale A^b, measured in units of g s a'N, is set by the vacuum expectation value 
of the dilaton $ . 

Numerically we are not doing anything new; the choice of the right endpoint of integration 
is always arbitrary, decided to best fit the physical problem at hand. Using a generalized 
shooting technique and integrating up to the point where Va& = VNonAb ± 0.00001, we find 

We thank Sean Hartnoll and Stathis Tompaidis for valuable dicussions and input regarding this para- 
graph. 



15 



Figure 2: VS(r) for abelian and non-abelian theories, dashed line and solid lines. 



an improved value for the WKB estimate, K = 4.291 (l/V g s oc'N) . This is a numerically 
stable eigenvalue meaning that small changes in the initital guess or pushing the endpoint 
of integration further to the right do not afect the value obtained. But this eigenstate is not 
normalizable, to obtain a normalizable state we have to impose the regularization procedure 
proposed above. This will be explained in detail in the next section. 

It is worth emphasizing that this model produces a discrete spectrum even without any 
kind of regularization. 

3.1 Understanding The Regularization Procedure 

Above we have proposed that the correct way of computing in this model is to do a compu- 
tation with the non-singular solution and, for large values of the radial coordinate, substract 
the result obtained with the singular background. This, we proposed is calculating in the 
dual N=l SYM theory. 

Let us get a better intuition of this sort of "regularization procedure" . In physical terms, 
this procedure is easy to understand and is just instructing us to do our computations only 
in the region that is of interest to N = 1 SYM. 

Indeed, since the non-abelian (non-singular) solution (that captures the IR effects of the 
dual field theory) asymptotes to the abelian (singular) solution (that is dual to a higher 
dimensional field theory), what we are basically doing when explicitely computing is substract 
the result obtained with the non-singular background minus those obtained at large values 
of the radial coordinate. Basically, this boils down to computing only in the region that is 
dual to N = 1 SYM (with KK impurities as explained above). This is not very different from 
the type of regularizations done, for example in the computation of Wilson loops [TT)].|17j. 
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where the infinite mass of the non-dynamical quark was substracted, or, what is the same, 
the mass of an infinite string not feeling the effects of the background (to which the real 
string asymptotes) is substracted. This sort of regularization was also used in [20] ■ In that 
paper, even when a hard cut-off was imposed for numerical convenience, one should think 
about it as the fact that the computation was done in the region of interest, where the probe 
brane that adds flavor to the quenched version of N=l SQCD is different from the probe 
brane in the singular background. More recently this sort of regularization was used in |24| . 
even when the model used in that paper is different from ours, we believe their regularization 
can be understood in the lines we wrote above. 

Some readers might object the following: if one computes in this way, everything will give 
a finite result, so in this case, all functions will be normalizable. This questioning is valid, 
so let us try to answer it; for this it will be convenient to resort on an example where an 
exact solution is known. 

Hence, it is instructive to analyze the solution studied in the paper [TT?] . Indeed, in that 
paper, the authors realized that a fluctuation given by 

%„ = ^=0, 5F 3 .F 3 = 0, 5F 3 = ndA 2 (3.37) 

with no restrictions to the functional form of the two form, solves the eqs of motion, that 
are written in the Appendix B (JB.lj) - (JB.8|) . This could be a massless glueball, but when 
computing the norm 

J d 10 x^\5F 3 \ 2 , 

seems to diverge, thus ruling it out as an state in the strong coupling theory. 

If we apply our criteria to this case, one might worry that the norm computed above will 
give a finite result, thus leading to a massless glueball that one does not expect in this theory 
(contrary to the KS case ^21)- If these quantities give a finite result, this will imply that an 
effect not expected (a massless glueball) shows up. 

So, to understand this, let us do the computation for the norm, and apply our regulariza- 
tion procedure 

\\5F 3 \\ oc / dre 2h+2(t> (3.38) 



the regularization proposed above, indicates that we do a computation like this 

\\5F 3 \\ = I" dre 2 * +2h \ non _ smgular + f°° dre 2 ^ 2h \ non _ smgular - [°° dre 2 *+ 2h \ singular (3.39) 

Where ro is a value of the radial coordinate where the functions e 2h+2< ^ computed in the 
singular and non-singular backgrounds are very similar to some degree of precision that 
is arbitrarily fixed. Since both integrands have the same asymptotics, they should equally 
diverge at large values of the radial coordinate. Indeed, both integrals diverge at leading order 
in the same way, but contrary to what one might expect, the integrals differ in a divergent 
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quantity (and many convergent terms), thus, the computation in ()3.39|) is divergent and the 
configuration in ()3.37j) is not a good state of our theory 

It is important to observe that many of the test that the solution has passed (see the 
review articles (23] ), still work with this regularization. 

We would like to stress that even what was done in ^3] was not technically correct (as 
we mentioned, they seem to have used the wrong fluctuated eqs), the hard cut-off that they 
introduced is doing the same job that the regularization that we proposed here. Nevertheless, 
we have to make clear some important differences with ^1] (apart from the fact that we use 
different eqs.). In [14] the authors did not find a discrete spectrum before imposing their 
hard cut-off, while we found one before our regularization, notice that our potential VS(r) 
does 'confine' wavef unctions. We need to appeal to the regularization, only to satisfy the 
normalizability condition of our discrete states (and this is because we are being conservative 
and adopting a curved space measure for our normalizations). So, the hard cut-off regulation 
is quite different from what we have done here. If one introduces a cut-off, together with 
some boundary conditions, all solutions to the Schroedinger eqs will be normalizable. In our 
case, things are more subtle, as we explained above. 

4 Glueballs from type IIA perspective 

This sections uses the same methods developed in Section 3 (See Appendix B.2 for all the 
details) for the type IIA background explained in Section 2. We will not carry out a full 
numerical analysis like in Section 3, but we will leave the system set for this more complicated 
(fully numerical) problem. 

Let us study glueballs for the case of wrapped branes in type IIA string theory. As 
explained above, in this case, the relevant background consists in N D6 branes wrapping a 
three cycle inside a CY3 fold. So, to summarize the things clearly, let us write the dilaton 
and the metric of the IIA solution in Einstein frame, 



The field strength F 2 reads, 
F 2 = k\r)drA(d^+cosOd<p-cosOd0)-(k{r)+l) sin 6d6Adp+(k(r) -1) sm9d0~Adip, (4.41) 




4e ^ = f(r) 2 + c(r) 2 (l+g 3 (r)) 2 . 



(4.40) 



with 



k(r) 



P-c 2 (i-gl) 

/2 + c 2(l + ^)2- 
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In the following, we will work with this IIA set-up and because of the many similarities 
with the IIB model studied in the previous section, we will use the same approach. It is 
interesting to mention that if we can find glueballs in IIA, they should have an expression 
in M theory purely in terms of a metric fluctuation. Let us start by finding the dynamics of 
the fluctuations. The Lagrangian of IIA is 5 

L = y/g(R- \{d^f - \e^Fl - ^Hj - le^ 2 F 4 2 ) + \b 2 A F 4 A F 4 , 

F 4 = dC 3 - H 3 A A u F 2 = dA 1 , H 3 = dB 2 , (4.42) 

now, let us focus on the configurations that are of our interest, that is those where only the 
fields 0, g^u, are turned on. The eqs of motion in this case are 

R,u = \d^d v <\> + l -e^(F^F: - ^„F 2 ), 

R = Imf + ± e 3/2*p2 V 2 0= 3 3 /20F2 

2 V YJ 16 ' Y 8 

d,(Vge 3 ^F^) = 0, d { ,F up] = (4.43) 

Now, let us study the variations of these eqs. Under a fluctuation in all the relevant fields 

(p = (p + e5(p, = + eh^ F^ = F^ + eSF^. (4.44) 

Again, we will propose a particular form for the metric fluctuation 

hij = ^rhg i:j , = x, t, r) h a/3 = h {aP) + -g a/3 . (4.45) 

The same comments that we made regarding this decomposition, before eq. (jrS.25)) are also 
pertinent here. We have denoted the trace of the internal part of the metric as h = h a a . We 
also propose an armonic expansion for the fields of a form similar to ()3.26|) . 

h(x,y) = Z I h I (x)Y I (y), h (af3) = Ejb^^Y 1 (y), 5<f>(x,y) = ^f 1 (x)Y 1 (y) (4.46) 

Keeping only the s-wave fluctuations as before, we find after many computations (that are 
carefully spelled out in Appendix B), that like in IIB case, the many equations nicely combine 
and is sufficient to solve just one equation that reads, 

V 2 h J + dg^d^drh 1 + — — (31F 2 + 57(g a PF a *F* - —g rr F r *F r *)) h 1 = (4.47) 

8 v H 57 ' 

As we have done in the type IIB section, the numerics in this case can be studied. We will 

not do this here and we just want to point to the fact that even when this has to be done in 

a completely numerical way (since the metric functions are only numerically known), there 

is a nice feature of this IIA solution. The fact that the dilaton does not diverge, makes us 

believe that the regularization will not be necessary. This is left for future work. The point 

of this section was just to call the reader's attention to this set of IIA models and show the 

analogy with the IIB treatment. 

5 notice that we take g s = 1 in this section 
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5 Summary, Conclusions and Future Directions 



Let us start by summarizing what we have done in this paper. First, we presented two 
Supergravity backgrounds (one in IIB and another in IIA) that are argued to be dual to 
N — 1 SYM at low energies and we analyzed the spectrum in detail. Then, we initiated the 
study of glueball-like excitations in the strong coupling field theory as fluctuations of the 
Supergravity fields. 

We presented the equations to study the spectrum of glueballs and their excitations and 
analyzed them numerically. The type IIB case is analyzed in detail and we proposed a 
regularization procedure that might be useful in other computations involving these wrapped 
branes set-ups. Two appendixes present our computations in full detail. 

We find that unlike some IIB backgrounds previously studied, in the D5 wrapped on S 2 
model not even the simplest scalar mode decouples from the rest of the fluctuations. Indeed, 
as we have shown, assuming only fluctuations of the dilaton leads to inconsistent equations. 
Therefore, the glueball ++ in the IIB model we studied, is not dual to the dilaton, but to a 
mixture of dilaton and trace of the internal part of the metric. This goes in the same line as 
the papers [Sj, where the glueballs turned out to be mixings between different Supergravity 
modes. This mixing might persist for higher spin modes. The presence of a non-constant 
dilaton background seems to be the reason for the mixing of the fluctuations and this also 
appears in the model studied in [12! ■ Another important point is that the potential found 
produces a discrete spectrum and a mass gap without any sort of cut-off, which seems 
to indicate that it is indeed capturing the physics of a confinig theory. As expected in 
a background with a linear dilaton, the states are not normalizable. Given that the UV 
completion of this theory is a little string theory it is not a surprise that a regularization 
(or substraction) procedure is needed. We present a proposal for this regularization which 
amounts to substracting the unwanted contribution of the UV regime. 

In the analysis of the type IIA background, we find that the scalar mode does not decouple 
from the rest of the fluctuations, indicating that, indeed, the non-constant dilaton in the 
background plays a role in producing this mixing. We do not perform the numerical analysis 
of the IIA background since the point of the paper is more to show a way of proceeding in 
these wrapped brane set-ups, that we believe is not exploited in the previous literature. 

Let us now discuss some future directions to follow and some work that should be inter- 
esting to do in detail, not discussed in this paper. 

First of all, as we mentioned in section 2, these models are contamined by the so called KK 
modes and we do not distinguish here if the glueballs we are obtaining are actually "made 
out" of KK modes composites (in which case they are not proper "glueballs" but hybrids 
composed of gluon, gluino and KK state) 

A good technical way to distinguish is to repeat the computation we are doing in Section 
3, but for the case of the "dipole deformed" field theory (see [13] for all details). Indeed, 
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the idea in the paper 03] is to make a SL(3, R) deformation of the supergravity solution, 
that reflects on the dual field theory side on particular deformation that affects only the 
dynamics of the KK modes in the spectrum. Hence, the comparison of the eqs (|3.27|) and 
those in the appendix with those obtained in using the deformed metric can illuminate on 
what type of composition our 'glueball' has, if only glue and gluino, or if it is composed out 
of KK modes. Same could be repeated with deformations of the G2 holonomy manifolds and 
the many examples already existing in the literature. 

It would be nice to check how our regularizing procedure works with the "flavor branes" 
introduced in j20j when doing a dual to quenched SQCD. Indeed, there the idea was precisely 
the same, taking away from the computation the effects of the unwanted UV region. This 
can be understood by looking at the plots in figures 2 and 3 of [20J. 

It should also be of interest, to study the glueballs in the type IIA model discussed in 
Section 4.1 of the paper jUj. This solution is basically the same as the one we discussed 
in the IIB section, after some dualities. So, the interest of studying glueballs in this case is 
obviously to see if the same spectrum is obtained, analyze differences among eqs of motion, 
etc. 

Besides, it might have some interest to study the IIA case in more detail, not only nu- 
merically, as we pointed out above, but also from a G2 holonomy perspective. Our dilaton- 
metric-gauge field fluctuations must combine in some way in a pure metric fluctuation in 
eleven dimensions. It should be nice to see how this works. 

Other models where it might have be interesting to apply our techniques (mainly the sort 
of manipulations explained in the Appendix B) is in models of non-supersymmetric duality. 
There is indeed one very clear model, that was studied in detail in the papers fjBj . 

One might also think about studying the 'fermionic counterpart' of what we have done in 
this paper, by fluctuating the fermionic fields around the bosonic background. 

On other respect, the comparison with the results from Lattice SYM should be done. 
There are some of these results in but we believe that the topic will evolve to allow better 
understanding and the contribution of this paper might be useful in the comparison with 
lattice results. See for example [47)|. Regarding this point, it should be interesting to study the 
spectrum of fermionic fluctuations (with fermionic fields vanishing in the background). This 
does not seem to have been exploited in the AdS / CFT literature, while other methods based 
on Veneziano-Yankielowicz and extensions, seem to give nice results [46J. This situation 
might clearly improve with some study. 

The interest of studying glueballs goes beyond the simple fact of getting a discrete spec- 
trum (that is by itself of enough interest). Indeed, glueballs play an important role in some 
recent advances regarding the study of Deep Inelastic and other types of Scattering using 
AdS / CFT techniques |UJ . The knowledge of glueballs masses and profiles in different models 
might help to extend the results in papers like [UJ to other 'more realistic' models. 
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A Appendix: Some Geometrical identities 

In the following, we list some geometrical identities that were used in the derivation of the 
fluctuated eqs flEHJ-jEEJ) 

6R = V^Vuh^ - V 2 K - R^ v 7 



B Appendix: Derivation of the eqs in the IIB and IIA 

cases 

In this Appendix we fill in all the details that were left out in the computations that lead to 
eqs. f|3.27|) and ()4.47|) in sections three and four. 

B.l Glueballs with the D5 branes solution 

Let us start with the type IIB solution. To study glueballs, we need the variation of the eqs 
of motion. In the type IIB case, for the solution of D5 branes wrapping S 2 , the Einstein eqs 
for the metric, dilaton and three form read, 



5( V 2 0) = V 2 64> - g ab 6T p ab d P 4> - h ab V a V h <P 




(A.l) 




(B.l) 
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by contracting we get the Ricci scalar eq, 



R=\^f+ 9 -^-Fl (B.2) 



The dilaton, Maxwell and Bianchi eqs. are, 



Fl dpL/ge+Fw) = 0, d^F K1/p] = (B.3) 



v 12 

Now, let us study the fluctuations of these eqs. above. Let us assume that the background 
fields vary according to 

SV«/ — ► g^ + th^, (f) -> + e<50, F^ p ^ F^ p + eSF^p (B.4) 

So, keeping only linear order in the parameter e, we get eqs for the fluctuated fields that 
read, for variation in the Ricci tensor eq. (|3.21|) . 

^[V a V,K + V a V„/£ - V„V^ - V 2 V] = \[d^dM + <9 A 



■r— [^F^F;* + 6F^F:* + F M **5F** - 2/i ab F^F; b ] 
-^[^(2F 3 5F 3 - 3h°»F a ^F b ** + 50F 3 2 ) + fc^Ff] (B.5) 
and for the Ricci scalar eq. (jB.2|) 

<IJ) [2F 3 5F 3 - 3h kl F k ^Fr + 50F 3 2 ], (B.6) 
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For the fluctuated dilaton eq. we have 



V 2 ^-/^V/A0-^^ = ( B - 7 ) 

and for the fluctuation of the Maxwell eq. and Bianchi identity, 



d^y/ge^d-hP + <50)F^ Q + 5F^ a - h ac F» v + h uc F? a - h^F™)} = 0, dSF = (B.8) 



1 

-/ 

•2 

Notice that the Ricci scalar eq. (|B.6|) can be obtained by contracting the Ricci tensor eq. 
(IB.5|) and substracting WR^, so, in the following we will work with a contracted version 
of (jB.5|) . In the derivation of these eqs, we have used the geometrical identities reviewed in 
Appendix A. 
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Now, let us assume a fluctuation for the metric of the form 



Mj ^9ijh a , (i, j x,t,T) 



Kp = h {a/3 ) + -g a ph a a , (a/3) = (0, 6, (p, <p, ip) 



(B.9) 



Where A is a constant and we have performed a shift in the h a p, this shift is proportional to 
the trace of the internal part of the metric /i". Note that the trace of the metric fluctuation 
over all the space is, htt = (A + l)h^. Now, let us study the form of the fluctuation of the 
eq. that is obtained by contracting (jB.5|) with the background metric g flu 5R llu , 



V„V^ - V 2 K = d r 8<j)d r (j) + 9 -^-m - Ig^h^F* + ZhTF^F? + 25F 3 ■ F 3 ] (B.10) 
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and using (|B.9|) . 
1 



AVX - (A + 1)(VJ + V>S - g"-d r 5<pd r 1 



24 



( A + 1 



+-(\ g "-F rm F? + g ab F a ^F b **)h a a + 2F 3 5F 3 + 3^^^ 
5 



2 
0. 



h«)F? 



:b.id 



For the dilaton eq. ()B.7jl we will have 



V 2 64>+(^±^)g"'d r( pd r h a a - 9se 



10 



12 



A. 
— / 

5 



+ '-hZ)F 2 + 2F 3 5F 3 - 



'h a a (\g rr F r ^F r ** + g ab F a ^F**) - Ah^F^F, 



(B.12) 



We have used the eq. of motion for the dilaton in the background eq. (jH.2Hj) above. The next 
step is to introduce an expansion in harmonics for each of the fields 

h%(x,y)=E I U I (x)Y I (y), h {aP) = E I b I (x)Y{ ap) (y), 6</>(x,y) = Zjf 1 (x)Y 1 \y) (B.13) 

In order to satisfy the Bianchi identity we write the fluctuation 5F 3 in the form 



5F„ 



d\JA v 



;b.i4) 



fivp [m v p] 

We will show that it is possible to find a fluctuation 5F 3 orthogonal to the background F 3 , 
i.e 5F 3 ■ F 3 = 0, that satisfies Maxwell's equation. 

Let us write the different componentes of Maxwell's equation for a general fluctuation 
5F 3 = d^dAyp] without demanding yet orthogonality with the backgorund F 3 . We get 

•u, a = a,b (angular) 



drlVge^d-h^ + 5<P)F rab + d [r 5A ab] - h bc F r c a + h ac F r c b - h rc Ff)\ 

+de[Vge 4 '{{-h^ + 5<f))F eab + 3 [e 5A ab] - h bc F e c a + h ac F 6b - h 6c Ff)] + d x [^(d [x A ab] )] = 

(B.15) 
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•v, a = r,b (r, angular) 



- h bc F e c r + h rc F 6b - h ec F r c b )} + d x [^(d [x 5A rb] )] = 



(B.16) 



*v, a = yz 



d r [^(d [r SA yz] )} + doWge+^Art)] + d x [^(d [x SA yz] )} = 0. 



(B.17) 



*u, a = (r, z) 



d 9 yge*{d^A r ^)\ + d^y/ge^d^SA^)] = 0. 



(B.18) 



Now demand 5F 3 • F 3 — 0, thus SF rab = 5F 9ab = 0. From equations (|B.15|) and (jB.16|) 
above it is clear that a fluctuation of the form 

- d x [^{d [x A ab] )} = d r [y/ge*{{~h» + 50)F rab - h bc F r c a + h ac F Tb - h rc Ff)} 



+ doyge+^h* + 5<P)F dab - h bc F° a + h ac F tib - h 0c Ff)} 



-d x [^ge^ x 5A rb 



■2~ i: 
1 



doWge%-h* + 5cj))F» rb - h bc F° r + h rc F° b - ff c F r c b )\ 



(B.19) 



with the other components given by ()B.17|) and (|B.18|) will satisfy Maxwell's equation, the 
Bianchi identity and is such that SF 3 ■ F 3 = 0. 

Therefore, we have only eqs. ()B.11|) and f)B.12|) . that keeping only the S-wave in the 
expansion in harmonics (jB.13)) take the form, 



and 



4A 

-(— + l)V 2 Mx) - d r f(x)d r $ - 
5 

+l(\g rr F r »F? + g ab F a ^Fnil(x) 
5 



,3A 



g s e v 

24 



(/(*) - ^±^n(s))i£ 







V 2 f(x) + (— + -)d r Tl(x)d r § - —e* (f(x) + -Tl(x))F 



10 



12 



A, 



(g ab F a ^F b ** + Xg rr F w F**)U(x) 







(B.20) 



(B.21) 



Equations ()B.20j) and (|B.21|) look suggestively similar so we will first check if there is any 
value of A for which they are the same. Indeed, for the two equations to be equal we need 
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to satisfy 



2( ^ + i)n(x) = /(*), 

2f(x) = Q + -)U(x), 

-(f(x) - { -^-U(x)) = f(x) + ^U(x) (B.22) 
It is easy to check that A = — 1| does the job. The eq. we need to solve is, 

a r5 2°. • 3 35 n 

V 2 f(x) + 2g rr d r f(x)d r <P+ ^\-Fl + — -(ff^W - -g rr F^F r **) ]f(x) = (B.23) 
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This is precisely the eq. f!3.27|) we wanted to obtain. 
B.2 Glueballs with the D6 branes solution 

The relevant background consists in N D6 branes wrapping a three cycle inside a CY3 fold. 
So, the metric of our IIA solution in Einstein frame and the dilaton where given in ()4.4()j) . 
and the Maxwell Field strength was, 

F 2 = k\r)drA(d^+cos6dip-cos6dip)-(k(r) + l)sm6d6Adip+(k(r)-l)sm6d6Adfi, (B.24) 

with 

k(r) = f'-^-ti) 

Because of the many similarities with the IIB model studied in previous sections, we will use 
the same approach. Let us first study fluctuations. The Lagrangian of IIA is 

L = MR - \d<P) 2 - -e 3 ^F| - -e~*H 2 - -e^F 2 
vy V 2 y 1 4 2 12 3 48 4 

+ ^B 2 A F 4 A F 4 , 

F 4 = dC 3 -H 3 A A u F 2 = dA u H 3 = dB 2 . (B.25) 

now, let us focus on the configurations that are of our interest, that is those where only the 
fields 0, g^, A^ are turned on. The eqs of motion in this case are 6 

^=1(50)2 + ^3/2^2 V 2 0= 3 3 /2 ^ 2 

2 V rj 16 Y 8 
dJygfWF^)=0, d [fl F up] =0 (B.26) 



D Like in the main text of the paper, we take g s = 1 in this Appendix 
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Now, let us study the variations of these eqs. Under a fluctuation in all the relevant fields 

= + e<50, g^v = g^v + e V = F^ + eSF^ (B.27) 
we have to first order in the fluctuation parameter e, 

SR^ = \{d^dJ4> + d v <t>dM) + ^^(^(F^F; - j^g^F 2 ) + F^5F* U + 

F U JF; - h ab F, a F ub - \g,v(F5F - h af3 F aif F;) - ^K V F 2 ), 

<5F = g^d^M - —^09,0 + -e 3 /^(-50F 2 + 2F5F - 2/^F^F;), 

5(V 2 0) = -e 3 /^(-<50F 2 + 2F5F - 2h^F^F*), 
8 2 

/ , F^ v \ 

^(v^ e3/ ^[( 35 + + - " ^"l) = ( B - 28 ) 

So, using the geometrical variations given in the Appendix A and putting all together, we 
end up with three eqs, one for the variation of the dilaton, one for the Ricci scalar and the 
Maxwell eq. One can see that the eq for the variation of the Ricci tensor is included in the 
Ricci scalar variation. 

So, we have 

V M V,/^ - V 2 /i£ = g^d^M + -g-f^^ 2 + 3F5F + ^F^F; - -^F 2 ) (B.29) 

a rr Q p 3/20 o 

V 2 50-^V M V,0-^-(2VX-V r ^)9 r 0= (-5<pF 2 + 2F5F-2h^F^F:) (B.30) 

^(v^e 3/2 ^[(350 + /i*)"2- + - K plM - KF au ]) = (B.31) 
Now, let us propose a fluctuation of the metric h^ u of the form 

hij = ^hgij, = x, t, r) h af3 = h {af3) + ^g af3 (B.32) 

Here, we denote the trace of the internal part of the metric as h = h a a . So, let us study the 
fluctuated eqs, we will have for the Ricci scalar, 

^5 5 ' 8^4 

+^F^F; + \g rr F r *F;)) - gW\V a V p h + F^Ff) = (B.33) 

and for the dilaton, after the eq. of motion for the background dilaton field has been used 
we will have, 

V 2 50 + (—±-)g rr d r <f>d r h - — {\^4> + ^h}F 2 + 6F5F - -h(g a ?F a *F; + \g rr F rif F*) 
—g^F^Ff) = (B.34) 

27 



and the Maxwell eq. 

^(v^e 3/2 1(350 + h)— + 5F^ - -g» r F? - g™hF£ - g^F? + 

\g» r hF» r - g» a hF: + g^Ffi) = (B.35) 

Now, let us follow an analysis very similar to the one we have done for the Type IIB case. 
First, we will propose an expansion of the form (|B.13j) . 

h{x,y)=X I h I (x)Y I {y), h (af3) = (x)Y( ap) (y), 84>(x,y) = ^f 1 (x)Y J (y) (B.36) 

Then, let us impose that both eqs are indeed the same. For this to happen, we will also need 
that the fluctuation of the Maxwell field is orthogonal to the field itself, that is 

8F^ = d^8A^ F 2 5F 2 = 0. 

(the first part is to automatically solve the Bianchi identity) and we will also concentrate on 
the s-wave, that is all modes with higher harmonics in the spheres will not be considered, 
same for g( a P> that contains higher harmonics. As before, the armonic decomposition is such 
that V^/i is composed of higher harmonics. Then dividing eq (|B.34|) by 6, we have that the 
following equalities have to be satisfied 

f =6(— g— )/i =(— )h , 3f = (— —)h, (B.37) 

for both eqs will be equal. Indeed, we can see that they are solved by 

A = -^, 95/ = 2/i, 
and that both eqs. (jB.33|) and (jB.34|) actually read, 

3/20 yi 

VV + §g TT d r $d r ti + — (31F 2 + 57{g a ?F a *F; - —g rr F r *F r *)) h 1 = (B.38) 
Regarding the Maxwell eq, we can make an argument similar to the one in the IIB case . 
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